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Abstract 

We analyze the LHC phenomenology of ASUSY — a version of NMSSM with a largish 
SH1H2 coupling. The scalar spectrum of the model contains a 200 — 300 GeV Higgs boson 
h with Standard-Model like properties, and heavy CP-even and CP-odd Higgs bosons H 
and A with masses in 500 — 800 GeV range. We study the discovery potential of H and A 
in the decay chains H ^ hh ^ AV — > 2/6j and A Zh Z2V '214: j. The dominant 
backgrounds are the diffuse Z6j and Z4j productions, which can be suppressed by demanding 
reconstruction of F's and h^s in intermediate states. The excess of signal events allows for a 
discovery of both H and A with over 5a significance for 100 fb~^ of integrated luminosity. 



1 Introduction 



The Naturalness problem of the Electroweak scale amounts to explaining the relative lightness of 
the Higgs boson compared to the Ultra- Violet (UV) cutoff of the theory. It is important to stay 
focused on this problem, since it provides the best hope to see new physics at the LHC. The scale 
of this new physics crucially depends on the mass of the Higgs boson: the lighter the Higgs boson, 
the lower this "naturalness cutoff' is expected to be. 

To make a quantitative estimate, one can use the indirect information contained in the Elec- 
tro Weak Precision Tests (EWPT). The Standard Model (SM) successfully accounts for the EWPT 
results for the Higgs boson mass rrih = 76^24 GeV [T] at 65% CL. The standard interpretation of 
this fact is that rrih should be quite close to its LEP2 direct lower bound of 114 GeV. To make 
such a light Higgs boson natural, new physics cutting off the top quark loop divergence should 
come in at or below the scale Anat — 400v^ GeV, if one allows finetuning of one part in A. 

The above standard interpretation of the EWPT has in it an implicit assumption that the new 
physics, while cutting off the top (and gauge boson) loops, does not itself contribute to the EWPT 
parameters T and S in a significant way. This assumption is, however, questionable, especially 
because the dependence of T and S on rrih is quite weak (logarithmic) in the SM. For example, a 
small breaking of the custodial symmetry in the new physics sector could be enough to generate an 
extra positive contribution to T making a much heavier Higgs boson consistent with the EWPT. 
This in turn allows to raise the naturalness cutoff of the theory by a non-negligible amount. This 
scenario, dubbed "Improved Naturalness", has been realized in simple explicit models [2] and 
leads to interesting modifications in the expected LHC phenomenology. 

Supersymmetry (SUSY), in many respects the most appealing way to UV-complete the SM, 
has its own specific problem with the Higgs boson mass. In supersymmetric extensions of the 
SM valid up to the GUT scale, such as the Minimal Supersymmetric SM (MSSM), the lightest 
Higgs boson h is generically predicted to be not much heavier than the Z. To achieve ruh > 114 
GeV requires a large radiative correction from the heavy stop, which introduces a few percent 
finetuning in the Z boson mass. This "SUSY Little Hierarchy Problem" led to several attempts in 
the literature to increase ruh by considering extensions of MSSM with extra sources for the Higgs 
quartic coupling. 

A representative example of these attempts is the Next-to- Minimal Supersymmetric SM (NMSSM), 
which contains a chiral singlet S interacting with the MSSM doublets Hi via a superpotential term 

XSH1H2. (1.1) 

The Yukawa coupling A contributes to the Higgs quartic coupling. Its maximal possible value 
is usually fixed by the requirement that it should stay perturbative up to the unification scale 
~ 10^^ GeV, in order not to disrupt the gauge coupling unification, which gives A^ax — 0.8 at the 
electroweak scale [3]. The resulting extra contribution to mh helps to bring it above the LEP2 
bound, somewhat reducing the needed finetuning |1] . However, the conclusion that rrih should not 
be much above 114 GeV remains unchanged. In fact this conclusion seems quite generic in most0 
extensions of MSSM which assume that all extra couplings stay perturbative up to the unification 
scale [H], [I!- 

notable exception is [S]. 
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On the other hand, if one abandons this assumption, rrih can be easily increased to a few 
hundreds of GeV. The most straightforward way to do this is based on adding the same term (11.11) 
to the superpotential, and taking A large. The scale Agtrong at which A becomes non-perturbative 
can be interpreted as the compositeness scale of (some of) the scalars, and the superpotential (11.11) 
can appear as an effective low-energy description of a confining SUSY gauge theory UV-completing 
the model above Astrong- 

The early studies of this idea [sj^ . which we call ASUSY to emphasize the role played by A, have 
shown that such UV-completions do exist and, moreover, can be made consistent with the gauge 
coupling unification. This positive existence proof is important in convincing us that ASUSY 
should be taken seriously. However, the details of the UV-completion are largely irrelevant for 
the TeV-scale phenomenology of the model, whose unique features are mostly determined by the 
presence of the large coupling A in the superpotential. 

Recently, the phenomenology of ASUSY, with the focus on the key issue of the EWPT, has 
been closely examined in [oj^ . and a very encouraging picture has emerged. Assuming that A 
remains perturbative up to about 10 TeV @, the lightest Higgs boson can be in 200 — 300 GeV 
range and yet consistent with the EWPT because of the extra positive contributions from the 
Higgs/Higgsino sector to the T parameter. These extra contributions are governed by the same 
coupling A as the Higgs mass, and thus do not require unnatural finetuning for cancellation. As 
a consequence of the increase in m^, superpartners such as stop can be in 500 — 1000 GeV range 
without finetuning (A = 5). This is the supersymmetric counterpart of the Improved Naturalness. 

Another attractive feature of ASUSY is the possibility of Higgsino Dark Matter. Due to strong 
mixing in the Higgsino sector induced by the same large A, the lightest Higgsino annihilation 
cross section is reduced compared to the MSSM case, and the observed Dark Matter abundance 
is reproduced in a large part of the parameter space [9]. 

All of the above makes ASUSY a well-motivated alternative to the conventional SUSY scenario. 
In this paper we would like to continue the study of ASUSY, focusing on its LHC phenomenology. 
Our purpose is twofold: we want to depict an LHC scenario which is impossible for conventional 
SUSY, but very natural for ASUSY, and we want to analyze how ASUSY could be observed at 
the LHC. We begin in Section [2] with a review of the model [9], stressing the differences with 
MSSM and NMSSM-like theories. In Section [3l we describe the early-stage LHC phenomenology 
of ASUSY, rather puzzling from the point of view of more standard SUSY scenarios. This puzzle 
would beg for an explanation and extra evidence, and we argue that it could come from observing 
the peculiar heavy scalars of the model. Our main result is a detailed study of the LHC discovery 
potential in the scalar sector of ASUSY at 100 fb~^ of integrated luminosity (Sections HI [5l). We 
summarize our conclusions in Section [6l 

^Originally called 'Supersymmetric Fat Higgs'. 

■^Scc also [To], where however the important Higgsino contributions to the EWPT have not been included. 
''This is necessary so that the EWPT can be analyzed in a fully perturbative way. 
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2 Review of ASUSY model 

The field content of ASUSY is tlie same as tliat of the NMSSM; that is, the only new field compared 
to the MSSM is a chiral singlet superfield S. The key feature of the model is the presence of the 
superpotential interaction f 1 1.11) with a large coupling A, which increases the mass of the lightest 
Higgs boson and improves naturalness of the theory, allowing for heavier superpartners. The 
maximal value of A is limited by the assumption that it stays perturbative up to about 10 TeV, so 
that the incalculable contribution to the EWPT from the cutoff can be neglected. In this paper, 
just as in [U], we take A = 2 at the electroweak scale. For this value of A the Landau pole is at 
about 50 TeV, which can be interpreted as the compositeness scale of (some of) the Higgs bosons 

2.1 Scalar sector 

The full ASUSY superpotential is 

W = fx{S)H,H2 + fiS), \ = fi'iS), 

while the scalar potential can be written in the form 

V = fil{S)\H,\^ + fil{S)\H2\^ - {fil{S)H,H2 + h.c.) + X^\H,H2\^ + V{S). (2.1) 

Here we neglected the gauge D-term contributions to the quartic term, which are small compared 
to the superpotential contribution for the chosen value of \U The mass parameters of the potential 
also include contributions from the soft SUSY-breaking Lagrangian. For simplicity, we assume 
CP invar iance of V and W. 

Many of the phenomenologically relevant properties of ASUSY can be characterized by the 
functions fij{S), fi{S), and M{S) = f"{S) evaluated at the Vacuum Expectation Value (VEV) 
s of the field S. These background values will be denoted below and M leaving their 

argument s understood. For example, the electroweak symmetry breaking (EWSB) is described 
by the equations 

tanp = — = — , 

AV = ^-/.?-/.^ (2.2) 
sm 2p 

where Vi^2 are the VEVs of the Higgs fields {v = {vf + v'^Y^'^ = 175 GeV). The mass of the charged 
Higgs bosons is 

The masses of the light neutral scalars can also be expressed via if their mixing with S can 
be neglected (we will comment about the validity of this approximation below). The mass of the 



^E.g.. the D-terms increase the mass of the hghtest Higgs boson of the model by 5 — 10 GeV compared to the 
expressions given below. 
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pseudoscalar A is then given by 




2 I \ 2 2 



sin 2/? 



The CP-even states hi have mass matrix 



m\ sin^/3 (A^^;^ - |m^) sin 2/? 

(A^t;^ — sin 2(3 m\ cos^/9 



The masses and compositions of the mass eigenstates /i, H are given by: 




(2.3) 



H = cos ahi + sin a /i2, h = — sin a /ii + cos a /i2 , 
_ m\cos2f3 + X 



(A^f^ — m?^±) sin 2/3 



It is convenient to parametrize the scalar sector of the model in terms of two parameters: tan (3 
and mii+. Their preferred range is: 



Here the bounds on tan (3 are suggested by the EWPT analysis [9] ; the lower bound on mH+ follows 
from requiring consistency with the constraint from b ^ s'-f without a destructive contribution 
from a stop-chargino loop [TT]; the upper bound on mH+ was derived in [H] from Naturalness 
considerations. The masses of neutral scalars in this range of parameters are given in Figs. [Tpi 
The key feature of the spectrum is that the lightest Higgs boson h is in the 200 — 300 GeV range, 
hence typically much heavier than in MSSM or NMSSM. Another notable feature (see Fig. [2]) is 
the fixed ordering of the spectrum: 



Throughout the paper we will assume that the singlet scalar S is significantly heavier than 
if^, H, A, and neglect the mixing between the neutral components of the Higgs doublets and the 
real and imaginary components of S. In fact S can be as heavy as 1 TeV or more, consistently 
with Naturalness [9]. In this case the no-mixing approximation works reasonably well. Analysis 
of concrete examples shows that the decrease in the masses of h, H, A due to their mixing with S 
does not typically exceed 5 — 10 %. At the same time the singlet admixture in h, H, A stays below 



2.2 Higgsino/Singlino sector 

In the fermion sector of ASUSY, we will assume that the electroweak gauginos are heavy and we 
will neglect their mixing with the Higgsinos. This is justifiable since the Higgsinos are mixed with 



1.5<tan/5<3, 
350 GeV < mH± < 700 GeV . 



(2.4) 



0.2-0.3. 
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Figure 1: The mass of the hghtest CP-even scalar h in the preferred region (j2.4p of the parameter space. 
The coupling A is fixed at A = 2. 
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Figure 2: The second CP-even scalar H and the CP-odd scalar A are always heavier than the charged 
scalars H^. This plot shows mass differences rriH — mH± (solid red lines, tan [3 = 1.5, 2, 2.5, 3 
from below up) and niA — fnfj± (dashed blue line) as a function of m[j± in the preferred 
region ()2.4p of the parameter space. The coupling A is fixed at A = 2. 
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the Singlino S by terms proportional to A, while mixing with the gauginos are controlled by the 
relatively small gauge terms [9]. In this case the charged Higgsino has mass fi. The neutral 
Higgsino mass matrix depends on fi and M; its expression is given in [9]. The lightest neutralino 
is always lighter than the chargino: 

m^o < m^± . (2-5) 
Stability of the potential gives an upper bound for the chargino mass: 

m^± < cos[3mH± , (2.6) 

which in turn implies that the lightest neutralino typically has a mass in 100 — 200 GeV range, 
so that it is the Lightest Supersymmetric Particle (LSP). Interestingly, this lightest neutralino 
can play the role of Cold Dark Matter. This is in contrast with the MSSM, where the pure 
Higgsino Dark Matter is disfavored since it (co) annihilates very efficiently and typically gives too 
low thermal abundance [12]. In ASUSY, the mixing between Higgsinos and the Singlino induced 
by A allow to reduce the annihilation cross section and get the correct thermal relic abundance. 



2.3 Other SUSY particles 

The masses of the top squarks and of the gluino affect the running of ^\ at the one- and two-loop 
level, respectively, and can thus be bounded from Naturalness considerations. For 20% finetuning 
(A = 5) and tan (3 as in (12. 4p these masses have to satisfy [9] 



mi < 600 - 800 GeV , (2.7) 
m-g < 1.2 - 1.6 TeV 

(looser bounds corresponding to smaller tan/3). For larger finetuning A these bounds increase by 
a factor a/A/5. 

The masses of the electroweak gauginos, sleptons and all the other squarks except for the stops, 
do not have significant Naturalness bounds. Thus it is relevant to consider the limit when these 
particles are well above a TeV. This limit is similar to the models with effective supersymmetry 
originally proposed as a way to address SUSY fiavor problems. 



3 ASUSY at the LHC 

We will now discuss the LHC phenomenology of ASUSY beginning from the easier signals of 
gluino, stop and the lightest Higgs and then continuing with a detailed analysis of the experimental 
signatures of the heavy scalars. In our discussion we will always assume A = 2 and tan j3 and mH± 
belonging to the preferred range (12. 4p . More specific Monte Carlo studies will be performed for a 
benchmark point 

mH± = 500 GeV, tan/5 = 2 , (3.1) 
corresponding to light neutral scalar masses of 

rrih = 250 GeV, niH = 555 GeV, = 615 GeV. 
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3.1 Gluino and stop 



The standard way to discover SUSY at the LHC is via pair-production of squarks and gluinos 

1221 [23]: 

PP 99, Q9, m ■ 

Since these sparticles are strongly interacting, the production cross section can be as large as a pb 
or more depending on the masses [I6]. The produced sparticles give rise to well-known cascade 
decays with lightest neutralinos in the final state, giving events with several jets, leptons and 
missing Et- 

The majority of available studies [221 [23] of this signal focus on the mSUGRA case, which 
gives degenerate squark spectra. While the same discovery strategy will apply also in the ASUSY 
case, the discovery is expected to be more difficult due to the fact that only stop squarks may 
be light enough to be produced. For a rough estimate we can use the existing study [14] of 
the LHC discovery potential in the case of effective supersymmetry [13], when only the 1st and 
2nd generation squarks are decoupled, while sbottom and stop masses are similar. Notice that in 
ASUSY the LSP is expected to be relatively light with respect to the gaugino and stop (see Section 
12. 2p . which helps the discovery. According to [H], in this favorable case 10 fb~^ of integrated 
luminosity should be enough for a discovery of SUSY in the entire range (12. 7p of stop and gluino 
masses suggested by Naturalness. 

The same signals with several jets, leptons and missing transverse energy can be used, in 
addition to SUSY discovery, to roughly estimate the sparticles masses. In fact, the total invariant 
mass of the visible particles in the final state cannot exceed the decaying sparticle mass. A concrete 
example is t — * llqxo] in this case one can set an upper bound on the stop mass 



where m^"'^ denotes the end-point of the invariant mass distribution of llq. 



3.2 Light Higgs 

The most peculiar property of the lightest Higgs boson h in ASUSY is its mass (see Figil]). This 
particle is always heavy (200 — 300 GeV), and this makes a relevant phenomenological difference 
with respect to the conventional supersymmetric models. To study other properties and analyze 
production and decay channels we need to know how the lightest Higgs interacts. 

According to the standard 2 Higgs Doublet Model (2HDM) result, the couplings of the h with 
the top and the weak gauge bosons are equal to the coupling of the SM Higgs boson times the 
following factors: 

cos d 

^htt = ^-^, ^hvv = sm{(3 - a) , 
sm p 

where from now on V means both Z and W. These are the only relevant couplings with SM 
particles. From Fig. [3]we see that the ^htt and ^hvv factors are very close to one (within 10%) in 
the whole parameter space of Eq. (12.41) . 

Since the h is SM-like, we expect the gluon fusion (GF) to be the dominant production process. 
For the moderate tan j3 of Eq. (12.41) the bottom loop is always negligible with respect to the top 
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350 400 450 500 550 600 650 700 350 400 450 500 550 600 650 700 



rriH' [GeV] ma' [GeV] 



Figure 3: (,htt (left) and ^hw (right), see Eq. ()3.2p . plotted in the range (\2A\\ for A = 2. 

loop. As discussed in Section [3.2. II below, correction due to the stop loop is also quite small. We 
conclude that the GF production cross section of h is always close to the GF production cross 
section of the SM Higgs boson of the same mass. 

On the other hand, decays of the h will also follow the same pattern as in the SM. In fact, 
the branching ratio is almost saturated by the decays into vectors. This allows to use existing 
SM Higgs boson studies to estimate the discovery potential of LHC. According to [221 ES], 5fb~ 
of integrated luminosity allows a 5a discovery in the "gold-plated" channel h ZZ —>■ . 
The mass of h will also be easily measurable thanks to good energy resolution of the final leptons: 
available studies [23] report that 30 fb~^ will be enough to measure rrih at 1 — 2 ppm level. 

3.2.1 Stop loop contribution to the GF production cross section 

The full (top+stop) LO amplitude for GF production of h is given by: 

AV^top = 9HuA^„,{n) + ^4(r,- ) , = £l • (3-2) 

i=l,2 ii « 

Here ii are stop mass eigenstates, and Qhi-i. are their couplings to the h. These depend on several 
parameters (stop soft masses mg, mu-, At as well as /i, a, (3) and are given by the same expressions 
as in the MSSM (see |2T], p. 24, 39). The and A^^ are, respectively, the stop and top loop 
amplitude (se [21], p. 92). At LO, the ratio of GF h production cross sections with and without 
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Figure 4: Kh plotted as a function of mq^u find the top trilinear SUSY breaking term Af. Other 
parameters are fixed as in Eq. (j3.4|) . The white areas in the lower left and lower right corners 
correspond to m^^ < 100 GeV (see Fig. |6]) and are therefore excluded by direct stop searches 

m. 



stop loop included is given by 



''"top+stop 
O"top 



Ah 

top+stop 



■^top 



(3.3) 



where by A^^^ we denote the first term in (13. 2p . To estimate the impact of the stop contribution, 
we evaluated for various plausible values of stop sector parameters. For example, in Fig. |4]we 
give a plot of Kh as a function of rhg u = {mij + mQ)/2 and At, with the other parameters fixed 
at 

tan /? = 2, mH+ = 500 GeV, tuq - mu = 100 GeV, /i = 200 GeV . (3.4) 

As we can see from Fig. HJ inclusion of stop loops changes the GF h production cross section by 
a small amount (less than 20%). We have checked that this result remains unchanged for different 
choices of the parameters in (13. 4p . The stop contribution is thus of the same order of magnitude 
as the NNLO QCD correction to the top loop GF process. 



3.3 What next? 

The early discoveries described in the previous two sections, if they indeed happen at the LHC, will 
be somewhat puzzling. Strongly-interacting cascade-decaying heavy particles will give a strong 
evidence for SUSY. At the same time, a SM-like Higgs boson with a 200 — 300 GeV mass will rule 
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out the MSSM, or other conventional SUSY scenarios. Indeed, in the MSSM, the hghtest Higgs 
boson mass has a theoretical upper bound of about 140 GeV lI, and there is no way to make the 
model compatible with the phenomenology of Sections 13.11 and 13.21 The same conclusion holds 
for all the extensions of the MSSM which keep couplings perturbative up to the unification scale. 

A natural way to resolve the puzzle is to allow couplings which become non-perturbative at a 
lower scale, and ASUSY is the simplest model which realizes this idea. At this point it will become 
crucial to make further tests of the model. With this in mind, below we will study the discovery 
reach of heavy scalars H and A and their mass measurement. This will be probably the simplest 
non-trivial measurement to perform. Interestingly, by measuring rriH and (and assuming that 
we have already measured nih), we can determine the main scalar sector parameters mH±, tsji (3 
and A (see Section [211]) . Knowing A, we can tell the scale at which compositeness/strong couphng 
sets in. 



4 The heavy CP-even scalar @ LHC 
4.1 Production 

The heavy CP-even Higgs boson H (see Eq. f l2.3l) ) has mass in the 500-800 GeV range (see Fig. [2]). 
Its couplings to fermions and weak gauge bosons are equal to the couplings of the SM Higgs boson 
times the following factors (see Fig. [5]): 

sin a cos a 

^mt = ^-5, ^Hbb = ^, ^Hvv = cos{fJ - a) . 

sm p cos p 

The Higgs-stop coupling gHi.ii is the same as in the MSSM and depends on fi, the a and (3 angles 
and on the top soft SUSY-breaking A-term At (see |21], pg. 24, 40 for explicit expressions). 

The relevant production processes are OF and vector boson fusion. Both top and stop loop 
contribute to the GF cross sectionj^ To estimate the relevance of stop contribution we studied 
the quantity K^, defined analogously to Kh from Sec. 13.21 In Fig. we plot and m^-^ as 
functions of friQ^u and At with the other parameters fixed as in fl3.4p . From this picture we can 
see that: (a) in most of the allowed area stop contribution enhances the cross section; (b) in the 
area corresponding to 

rui^ > 400 GeV (4.1) 

the stop loops correct the GF cross section by less than 20%. As in the light Higgs case, such a 
contribution is comparable to NNLO QCD correction to the top loop diagram, which is never taken 
into account in this paper. In what follows we will neglect the stop loop contribution to the H 
production. Thus, in most of the parameter space we will be underestimating the production cross 
section. In principle, as discussed in Section IXTl stop mass is likely to be determined/constrained 
from cascade decays. Therefore we would be able to check assumptions like (14. ip and improve the 
accuracy for GF H production cross section prediction. 

*5 Assuming mstop < 2 TeV, see [21]. 

^The bottom loop contribution can be safely neglected because £,Hbb is not large enough to make it comparable 
with the top loop for moderate tan /3 as in (|2.4p . 
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Figure 5: ^//(^(left) and ^nvv (right) plotted in the range <\2Ah for A = 2. 



We thus obtain the gluon fusion (GF) and vector boson fusion (VBF) production cross sections 
of the H by simply rescahng the NLO results for the SM Higgs boson of the same mass, generated 
by HIGLU [17] and vv2h [18] codes. For instance the GF result is given by: 

„GF ( TT\ £-2 GF 

The obtained production cross sections are shown in Fig. [71 With an order of 0.1 pb production 
cross section (GF being the dominant mode), the search for H looks feasible. 

4.2 Decays 

The visible decay width of the H is dominated by decays into hh, tt and VV pairs. We assume 
that the decay channel into a stop pair is closed, which happens in most of the parameter space 
(see Fig. [H]); such an assumption can be checked and corrected if necessary when the stop mass is 
measured. The Hhh coupling is proportional to and is given b}|§ 

guhh = —7= [sm(a + /5) - 3 sm(3a - p)] . 
l\j 2 

The total visible decay width is given in Fig. [H] and ranges between 5 and 25 GeV. The branching 
ratio for decays into hh and, for comparison, into ZZ pairs, is plotted in Fig. Because A is large. 



^The corresponding Lagrangian term is guhkHh? /2. 
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Figure 6: Left: mass of the lightest stop plotted as a function of rriQ^jj and Af. Right: The Kh ratio 
plotted as a function of ffiQ^jj and At- In both plots parameters are fixed as in Eq. (j3.4|) . 
The white areas in the lower left and lower right corners correspond to nif < 100 GeV and 
are therefore excluded by direct stop searches [T^ . 
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Figure 7: The NLO production cross section of the H via the gluon fusion (left) and the vector boson 
fusion (right) plotted in the range ()2.4p for A = 2. 
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Figure 8: Th omitting supersymmetric decays into Higgsino pairs, see Fig. [TUJ 



decay into hh pairs is a dominant decay mode whenever this channel is open, which happens in 
most of the parameter space. 

The H will also decay into Higgsino pairs. This decay width depends on the Higgsino sector 
parameters /i and M , see Sec. 12. 2[ Fig. [10] gives the decay width of H into Higgsinos for the 
benchmark point (13. ip and for /i, M within their ranges (determined by stability of the potential 
and Naturalness considerations [9]): it takes values between a few and 15 GeV. Below we will 
neglect the decay width into Higgsinos. This means that in a realistic situation all branching 
ratios and signal rates will have to be multiplied by a factor r/(r + T^^). 

4.3 Detection strategies 

Let us first discuss the lower left corner of the parameter space, where the H hh decay channel 
is closed (see Fig. IH]). In this region BR{H VV) becomes significant. We believe that H 
could be discovered in this region via H ZZ 41, vvll combined with H —>■ WW — > lulu. A 
rough estimate of the discovery reach can be obtained using results of the SM Higgs boson studies 
[221 [23] . and then taking into account that the width of our H scalar is significantly smaller that 
the width of the SM Higgs boson of the same mass. More precisely, the discovery significance can 
be estimated by rescaling the corresponding significances in the SM case with a factor 
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Figure 9: BR{H hh) (left) and BR{H ZZ) (right) in the preferred range (12. 4p of the parameter 
space. The H decay width into Higgsinos F^^ is neglected. For nonzero F^^, these branching 
ratios have to be multiplied by a factor F/(F + F^^), where F is the visible decay width 
plotted in Fig. El The H ^ hh decay mode is dominant except for the lower left corner of 
the parameter space where this decay channel is closed [mn < 27?t,/j). 




Figure 10: The H decay width into Higgsino pairs for A = 2, and tan/3 at the benchmark point 

(13. ip . and for fj, (chargino mass) and M within their ranges determined by stability of the 
potential and Naturalness considerations [9]. The gray area corresponds to mLSP < wz/2 
and is excluded. 
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where the quantities marked by SM refer to the SM Higgs boson of the same mass as the H. The 
factor \/Tsm/^h reflects the reduction of background events passing the event selection in the 
mass window ±const.T. This rescaling procedure gives a 5 — 60" significance with 100 fb~^ for the 
H discovery in H ^ VV when H ^ hh is closed □. 

In the remaining, larger region of the parameter space, BR[H —>■ VV) is too small for a 
convincing H discovery in the VV decay channel. In what follows, we will discuss how H could 
be discovered in that region using the decay mode H hh. The fact that this decay mode is 
dominant when open reflects a very basic property of ASUSY: the large value of A. 



4.4 Signal from H ^ hh 

For H ^ hh decay we cannot rely on existing SM studies. To perform a careful analysis, we 
will consider a benchmark point (13. ip . This point is generic rather than chosen for some special 
properties. The relevant particle parameters at this point take the following valued: 

ag^ = 150 fb, = 27 fb, 

mn = 555 GeV, nih = 250 GeV, 

Th = 21 GeV, Th = 3.8 GeV, (4.2) 

= 0.058, envv = 0-060, 

BR{H hh) = 0.76, BR{H ^ VV) = 0.2 . 

As discussed in Section 14.11 the H is mainly produced via gluon fusion; in the following we 
will consider only this channel. Once produced, most of the Hs will decay into hh and then 
into iV, resulting in agg^n^w = 110 fb. The final weak bosons can decay leptonically, but the 
branching fractions in this case are too small to allow more than one leptonic decay. Our choice 
for a quantitative study is the channel with one leptonic Z decay, with the remaining weak bosons 
decaying hadronicalljcj: 

gg ^ H ^ hh ^2Z2V ^ l-^reJ , a x BR = 2.67 fh. (4.3) 

To increase the signal cross section, we assumed that final state jets J are generic jets, i.e. 
J = j, b, c, where j is a usual gluon or light-quark jet. Flavor labels are not necessary, since we 
will not deal with flavor tagging issues at all. 

To produce a sample of signal events, we first used madgraph [21] to produce matrix-element- 
generated gg ^ H ^ VVZl'^l' events, and then we simulated the decay of the remaining weak 
bosons through the DECAY routine by F. Maltoni 



^The preceding discussion used the gluon fusion production mechanism. This result can presumably be improved 
using vector boson fusion, which is not normally used in the SM for this range of the Higgs mass, but becomes 
significant in ASUSY for low mjj+ (see Fig. [T]). 

^°The reported decay widths and branching ratios are calculated assuming zero decay widths into Higgsinos. See 



discussion in Section [ 

^^The alternative channel H — > WWW — > IfQJ benefits from a higher rate and could perhaps yield a higher 
statistical significance. Another promising channel is H ^ WWWW with several same-sign or opposite sign- 
different flavor leptons in the final state, which was recently used in a related study of non-SUSY H ^ hh decays 
[25j . We preferred channel (j4.3p to avoid discussing additional sources of missing energy among which there are 
particularly delicate detector effects (jet energy scale, finite cone size effects, calibration, etc.). 
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4.5 Backgrounds 

We scanned the long list of SM processes with l~^l~6J final state and used alpgen[29] or mad- 
GRAPH+DECAY to compute their cross sections for the total invariant mass near the H mas^. 
We found that only Z6J and ttZ processes are relevant, i.e. have cross section large enough to 
potentially compete with the signal. The details of this preliminary analysis can be found in 
Appendix \M 

We then proceeded with a more complete analysis of these two relevant backgrounds. Samples 
of {Z — > l~^l~)6j and {Z — > l'^l~)AiQQ {Q = c,b) events were generated with alpgen using the 
CTEqSl parton distribution functions (pdf). We used cuts 

ARjj>0.7, p^>20GeV, r/j < 2.5, (4.4) 
80 GeV < ma < 100 GeV, rji < 10. 

We also enforced the total invariant mass cut 

400 GeV < mtot,in^ < 2400 GeV, (4.5) 

covering by a large margin the region near the H mass. This allows us to properly introduce a jet 
spectrum smearing and take into account possible effects from high invariant mass tail£f|. With 
these cuts, our results for the cross section are reported in Table [TJ These results were obtained 
for the ALPGEN factorization and renormalization scale set at = m| + Pt z- motivation 
for choosing this scale is twofold. First, the Tevatron experiments [301 El] have confronted the 
observed rates of Z + N jets events with alpgen simulations for various fi'jp, finding fijp values 
not far from our choice as best fitting the observations. Second, our yields nearly the largest 
cross section we found trying out several possibilities available in alpgen. Thus we believe that 
the systematic uncertainty of background normalization is conservatively taken into account. 

The {Z l^l^)ti process, with subsequent 6J decay of the ti pair, was simulated with 
MADGRAPH+DECAY using the CTEqGlI pdf. We generated a sample using cuts (14. 4p and setting 
the renormalization and factorization scale at = m|; see Tabled] for the cross section estimate. 

4.6 Analysis 

The total background cross section, see Table [H, is much bigger than that of the signal, eq. (14. 3p . 
However, we expect signal events to have very specific structure due to the presence of intermediate 
resonances {h, W, Z). Typical background events are not expected to have such structure and can 
be rejected by imposing reconstruction cuts, i.e. requiring that the intermediate state resonances 
be reconstructed by final state jets and leptons. This is the general idea of the analysis described 
below. The main issue is whether the rejection efficiency will be enough to sufficiently suppress 
the backgrounds. 

^^Possible SUSY backgrounds like sparticle mediated diffuse hh production and l^l^QJ + LSP in gluino and 
squark decay have been estimated to be negligible. 

^•^At the same time, the lower invariant mass cut in ()4.5p was indispensable with our limited computer resources, 
since it improves greatly alpgen unweighting efficiency. 
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Process 


specific cuts 


a 


{Z i+i-)Qj 




1118(2) fb 


{Z l+r)bbAj 


p't > 10 GeV 


94(2) fb 


{Z l+l~)ccAj 


pir > 10 GeV 


92(1) fb 


{Z l+l~){tt 6J) 


rii < 2.5,pV > 10 GeV 


5.86(2) fb 



Table 1: Simulation of the relevant SM backgrounds for H l~^l~6J. Apart from the shown specific 
cuts, all pairs and jets fulfill Eq. ()4.4I) . Final state total invariant mass is between 400 
and 2400 GeV except for Ztt, which is produced without invariant mass restrictions. 

The details of the analysis (performed with root|26]) are as follows. First of all, our analysis is 
completely partonic, so neither showering nor jet reconstruction effects are taken into account. We 
also ignore flavor tagging and trigger issues, but our inclusive deflnition of jets and flnal selection 
cuts for leptons make these simpliflcations fully justifled. However, in order to make the analysis 
more realistic, we do introduce a smearing of energies of individual jets using the expression 

to generate the smearing coefficient, as prescribed in Table 9-1 of [28 



After smearing, we impose the kinematical cuts (14. 4p on the jets and slightly tighten the 
corresponding leptonic cuts 

ARij > 0.4, plj. > 10 GeV, rn < 2.5 . (4.7) 

Background and signal events not passing these cuts are removed from the samples. The signal 
events passing these cuts correspond to the 0.42(1) fb cross section (to be compared with Eq. 
(14.31) without any kinematical cuts). The background cross sections are reduced by these tighter 
cuts only by a small amount compared to the values reported in Table [H 
Finally, we impose the reconstruction cuts, proceeding as follow£f|. 

Rl. For each event we try to group the 6 final jets into 3 pairs so that the jets in each pair 
reconstruct a W or a Z. By this we mean that the invariant mass rriinv of each pair has to satisfy 
the requirement: 

Mv-6v< m,nv <Mv + 6v, 6v = 8 GeV, V e {W, Z} . 

In practice, the value of 6v cannot be taken too small because otherwise too many signal events 
will be rejected. The given value was motivated by the finite resolution of the W and the Z peaks 
which is determined by their natural widths as well as by the energy resolution of the detector as 
taken into account by the smearing procedure described above. 



^^This is also the smearing adopted in the ATLFAST++ detector simulator p7] . 

"'^^ Geometrical discrimination has been attempted too, but turned out not to be very helpful, since both signal 
and background result in a largely boosted system. 
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Process 


a 


H GJl^l' 


0.286(8) fb 


{Z l-l+)Q3 


0.15(1) fb 


(Z ^ l-l+)QQAj 


0.028(3) fb 


{Z ^l-l+){tt 6J) 


0.022(1) fb 



Table 2: Signal and background cross sections after imposing the reconstruction cuts. For ZQQAj the 
given value is the sum over Q = b,c. 



R2. If a grouping into jet pairs reconstructing a W or a Z each is found, we proceed to impose 
a further condition that two /I's be reconstructed by four jets from two of these three pairs, say 
pair 1 and 2, and by two jets of pair 3 and the two leptons. In this case the precise reconstruction 
cut that we used is 

mh-Sh< rUpairi+pairi < T^h + ^h, = 18 GcV, 

rrih- rripairs+i+i- <mh + —, 

where 'mpairi+pair2 and mpairs+i+i- are the invariant masses of the 4J and 2Jl^l~ final states. The 
value of 6h is again motivated by the natural width of h (with an additional spreading caused by 
the jet energy resolution). We also check that the gauge boson reconstructed by the jets of pair 3 
is a Z, while the two gauge bosons reconstructed by the jets of pairs 1 and 2 have the same type 
(both W or both Z). 

If no grouping of 6 jets into 3 pairs satisfying both Rl and R2 can be found (we go over 
all combinations), the event is rejected, otherwise it is retained. The retained events show the 
expected intermediate state resonance structure of the signal. 

We ran the reconstruction analysis on the signal sample and on each of the relevant background 
samples shown in Table [H The signal and background cross sections after the reconstruction cuts 
are given in Table [21 For each sample the number of events which passed the reconstruction 
cuts was large enough so that the statistical uncertainty in determining the rejection efficiency 
is reasonably smalL£|. In fact it is this statistical uncertainty (determined from the usual 
fluctuations of the number of events passing the reconstruction cuts) which underlies the errors 
for the cross sections quoted in Table [21 

Two basic conclusions are evident from Table [21 On the one hand, we see that for the chosen 
parameters 5v,h the reduction in the signal cross section from what we had after the kinematical 
cuts is reasonably small (from 0.42 to 0.28 fb). On the other hand, we see that the reconstruction 
cuts have huge effect on backgrounds, giving the rejection efficiency of about 10~^. The final 
background cross section is comparable to that of the signal, making the discovery possible. 

Finally, in Figure [TT] we show the distribution of the signal and the total background cross 
section versus the total invariant mass of the event. 

^^For example, we had 1059 events in the signal sample which passed all the cuts. 
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SoO 450 500 550 600 650 700 750 800 850 900 

Mi„JGeV] 

Figure 11: ASUSY at the benchmark point (13. ip . ()4.2p (black) and Standard Model (grey) expectation 
for the differential cross section da/dMinv{l'^l~QJ) after the kinematical and reconstruction 
cuts discussed in Section 14.61 

4.7 Discovery potential after 100 fb~^ 

From Figure [TT] we see that signal and backgrounds peak in the same invariant mass range. The 
discovery of H will thus come not from an observation of a new peak, but rather from an overall 
excess of events compared to the SM prediction, as well as from the enhanced prominence of the 
SM peak. 

For an integrated luminosity of 100 fb~^, the expected number of events passing all the cuts 
is 20 in the SM, and 49 in ASUSY at the benchmark point (13. ip . (14. 2p . giving SAa if one uses the 
significance estimator formula given in Eq. (A. 3) of [23]. Of course, once this global excess is found, 
it is worth to scan the invariant mass range to find where the excess is localized. Optimizing the 
range, much better discovery significance can be achieved. For instance, for 510 GeV< Mj„,„ < 590 
GeV we have 4 events in the SM, and 24 events in ASUSY, 6.860" away from the SM. When going 
beyond benchmark-point analysis (something we do not attempt in this paper), such localized 
excess can be used to determine niH- 

Our conclusion is that the ASUSY signal (14. 3 p is indeed observable at the LHC with 100 fb~^ 
of integrated luminosity. If observed, it can provide clean evidence for the heavy scalar H as well 
as for the H ^ hh dominant decay chain. 



19 



a(gg^A) [pb] 




350 400 450 500 550 600 650 700 
mn^ [GeV] 



Figure 12: Pseudoscalar Higgs boson production cross section plotted in the parameter space of Eq. 
m for A = 2. 

5 The CP-odd scalar @ LHC 
5.1 Production and decays 

The pseudoscalar Higgs boson A has mass in the same 500-800 GeV range as the heavy scalar H, 
but is always heavier than H (see Fig. [2]). Its couplings to the third generation SM fermions are 
given by [21] : 

gAtt = — cotp , QAbb = — tanfj. 

V V 

i > 

There is also an gAhz{A d ^h)Z^ coupling [21] 

2 cos t)w 

where g is the SU(2) gauge coupling. 

By CP-invariance AVV couplings vanish, therefore the only relevant production mechanism 
of A is gluon fusion. GF cross section is dominated by top loop. The bottom contribution is not 
significantly enhanced for moderate tan (3 and remains negligible. Stop loops are absent due to 
the combined effect of CP, which forbids Atiti and At2t2 couplings, and gauge invariance, which 
forbids gluon couplings to tit2 states. 

The production cross section has been evaluated at NLO with HIGLU [H] and is plotted in 
Fig.m 
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BR (A^hZ) 
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Figure 13: BR(74 — > Zh) plotted in the parameter space of Eq. (|2.4|) for A = 2, assuming negligible 
decay width into Higgsinos. BR{A 1-BR(A Zh). 

The total width of A ranges between 5 and 30 GeV and is dominated hj A tt and A — > hZ 
decays. Although the branching ratio of A tt is almost always dominant (see Fig. [T^ . we 
cannot exploit this channel. Indeed, [12] showed that for the mass values we are interested in, 
the ti SM background does not allow discovery of a scalar resonance decaying into ti. Therefore, 
we focus on A ^ hZ, whose BR is smaller, but still significant. Most of the produced /I's will 
decay into vectors, yielding (Jtot{gg ^ A ^ ZVV) ~ 100 fb over all the parameter space. Such 
a cross section will give too small event rate if more than one V is allowed to decay leptonically. 
Therefore we concentrate on the signature 

gg ^ A^hZ ^ VVZ AJl+r (signal). (5.1) 

For a detailed study we go to our benchmark point (13.11) . which gives the following numerical 
valueJ^: 

mA = 615GeV, = 11 GeV, (5.2) 

a X fii?(signal) = 6.9 fb. 

^^The quoted value of Ta does not include the width into Higgsino pairs depending on fi and M. The latter can 
be as large as 10 GeV, but in most of the parameter space is below 2 GeV. See the analogous discussion for the H 
in Section [421 
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Lnannel 


a 


A / ^ ^ T — 1— 7 \ A ^ 

[Z ^ /+/ )4J 


3.02(4) fb 


{Z l+l-)AJ 


7.006(4) pb 


{z in-)w2j 


176.0(8) fb 


Sum of neglected 


~ 90 fb 



Table 3: Cross sections of the signal and of the relevant SM backgrounds after the kinematical cuts 
()5.3p . Notice the reduction in the signal cross section compared to ()5.2p . 

5.2 Analysis and discovery potential 

The analysis is quite analogous to what was done in Section H] for the heavy CP-even Higgs. We 
will therefore be relatively brief. 

We generated a sample of signal events with madgraph+decay. 

We then computed cross sections for all SM processes with AJl~^l~ final state and we found 
that only the ZAj and ZW2j processes are relevant^. Table [3] contains details about the relevant 
and neglected backgrounds. 

Event samples of relevant backgrounds were generated for a more complete analysis. The ZAj 
process was simulated with alpgen using the CTEq5l pdf, setting fij. = m| +Px z enforcing 
cuts: 

500 GeV < Mi„„ < 750 GeV (5.3) 

> 20 GeV, > 10 GeV, 
ARj.^ij^a > 0.4, T]''' < 2.5 
80 GeV < mil < 100 GeV . 

The ZW2j process was simulated with madgraph using the cteqGlI pdf, setting = m| and 
imposing the same cuts. The resulting background cross sections are the ones given in Table [3l 
Finally, to model finite detector energy resolution, we apply energy smearing to the signal and 
background events using the smearing function 04.61) . 

From Table 3 we see that the background cross section in the relevant interval of rriinv is more 
than 3 orders of magnitude above the signal. To reduce the background, we proceed by imposing 
the reconstruction cuts. Namely, we require that the 4 final jets can be divided into 2 pairs 
reconstructing a vector boson each. Moreover, we require that these two vector bosons be of the 
same type. If they are both W, then we require that they reconstruct an h. If they are both Z, we 
require that out of the 3 final Z's (the two from jets and the one reconstructed by the leptons) we 
should find two reconstructing an h. Reconstruction parameters 6v,h, having the same meaning as 
in Section SSI are taken to be 

Sv = 8GeV, 4 = 18 GeV. 

^®In particular, backgrounds ZZ2j, VVV, h2j were analyzed and found negligible. 
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Figure 14: ASUSY (black) and SM (grey) expected differential cross section da / dMinv{4:jl~^ l~) for pro- 
cess ()5.ip at the benchmark point (|3.ip . (j5.2p . 

If the above requirements can be satisfied, the event is retained, otherwise it is rejected. 

The portion of the signal event sample which passes the reconstruction cuts amounts to 2.2 fb; 
cross section is reduced only by a small factor compared to the value after kinematical cuts given 
in Table [3l At the same time the total SM background cross section in 500 — 750 GeV invariant 
mass range drops after the reconstruction cuts by a factor of about 200, to 51.1 fb. 

Of interest are the differential cross sections of background and signal+background versus the 
invariant mass, plotted in Fig. [T31 We see that the signal distribution presents a clearly visible 
peak above the background. The discovery significance can be optimized choosing a range with 
largest S/\fB ratio. For example, assuming 100 fb~^ of integrated luminosity, in the 595 — 635 
GeV range we expect 816 events in the SM, and 989 events in ASUSY at the (13. ip . (15.21) benchmark 
point, which amounts to 6.I0" discovery significance. 

In summary, we showed that the CP-odd Higgs boson of ASUSY has a clear experimental 
signature in the 4j7^/~ channel, allowing for its discovery at the LHC with ~ 100 fb~^ of inte- 
grated luminosity. Moreover, the peaked shape of the signal distribution should allow background 
extraction from data and an easy mass measurement. Even though the A —>■ Zh decay mode is 
less distinctive of ASUSY than the H —>■ hh mode discussed in Section HI its signature is much 
simpler and cleaner, and it could be the easiest channel to pursue when looking for ASUSY signals. 

6 Overview and conclusions 

Soon the LHC will start directly probing energies well above the electroweak scale. If low-energy 
SUSY is the mechanism which stabilizes the gauge hierarchy, it should be discovered by the LHC 
experiments. However, the parameters and even the full field content of the fundamental La- 
grangian will be much more difficult to determine than the existence of supersymmetric particles. 
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In this respect the mass of the hghtest Higgs boson can be an interesting indicator. In the minimal 
version of supersymmetry, the MSSM, the Higgs boson is generally predicted to be light with a 
theoretical upper bound of around 140 GeV0. Thus the observation of an heavier Higgs boson will 
rule out the MSSM and all its extensions which keep couplings perturbative up to the unification 
scale. 

A well-motivated alternative which may realize the heavy h scenario is ASUSY [9]. The key 
feature of this model is the introduction of a chiral singlet superfield coupled to the MSSM Higgs 
doublets by a cubic superpotential term XSH1H2. Since the coupling A is taken largish, A ^ 2 at 
low energies, rrih is naturally in the range 200 — 300 GeV. 

If ASUSY scenario is realized in Nature, we expect early discovery of SUSY via gluino and 
stop decay cascades and the discovery of a SM-like lightest Higgs boson of mass around 200-300 
GeV. These experimental results would be at odds with the lore of SUSY phenomenology. They 
are, however, very natural in ASUSY. 

In this paper, we have investigated how to continue the experimental study of ASUSY if this 
puzzling picture emerges after 10 fb~^ of integrated luminosity will be available. The plan is to 
search for the heavy neutral scalars H and A, which are expected to be in the 500 — 800 GeV 
rangj^. More specifically, we have studied a benchmark point (13.11) where rriH = 555 GeV, 
rriA = 615 GeV, while the mass of the lightest Higgs boson is rrih = 250 GeV. The production 
cross sections of H and A are 177fb and 0.7pb, respectively, and are dominated by the gluon 
fusion. We have studied the discovery reach in the decay chains 



The dominant background for these signals are the diffuse ZAj and ZQj production. We suppress 
these backgrounds by demanding reconstruction of the vectors and the light Higgses in the inter- 
mediate state. The excess of signal events over the suppressed background allows for the heavy 
scalar and pseudoscalar discovery with over 5a significance for 100 fb^^ of integrated luminosity. 
Moreover, we found that A's background shape can be extracted from data, which should lead to 
a discovery with even less than 100 fb~^. Notice that the heavy scalar signature, relying on the 
Hhh coupling, is particularly distinctive of ASUSY. 

The masses of both H and A can be roughly determined by the position of the peaks in the 
total invariant mass distribution. Their masses (together with the early measurement of rrih) will 
allow to determine the 3 parameters of the Higgs sector: mj|/+ , tan /5, A (see Sectior l2?T]) . The 
determination of A is especially important since this coupling is the cornerstone of the ASUSY 
hypothesys; measuring A will also allow to evaluate the scale at which the non-perturbative physics 
sets in. 

It would be interesting to find which further checks of ASUSY could be performed by the 
LHC. The measurements of h GF production cross section could be one such check. Since h is 
SM-like, its production cross section in the interesting rrih range can be determined with a 10% 
accuracy using 100 fb~^ [221 123]. On the other hand, the theoretical uncertainty is around 20% 
and basically arises from the ignorance of the stop loop contribution (see Section [H^ . 

-"^^ Assuming mgtop ^ 2 TeV, see pT] . 

^'^Thc fourth heavy scalar, predominantly a singlet, can have a mass above a TeV and is unlikely to be observed. 



A^hZ 2VZ 4:jl+l-, 



(6.1) 
(6.2) 
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Process 


specific cuts 


a 


{Z ^ /-/+)6j 




34.2(2) fb 


{z ri+)bb4j 


T 

p't > 10 GeV 


4.22(2) fb 


{Z l-l+)ccAj 


pir > 10 GeV 


4.0(1) fb 


{Z l+l-){ti 6J) 


r]i < 2.5,pfr > 10 GeV 


5.86(2) fb 


Sum of neglected 


See Appendix "X 


^ 1 fb + ZVAj < 12 fb 



Table 4: Preliminary simulation of SM backgrounds for H I'^l^dJ used to select the two relevant 
backgrounds: Z6J and ttZ. Apart from the shown specific cuts, all pairs and jets fulfill 
Eq. ()4.4p . Final state total invariant mass is between 500 and 600 GeV except for {Z l+l-)tt, 
which is produced without invariant mass restrictions. 

A more non-trivial test of the theory could come from the direct observation of H~^. Whether 
this is feasible needs a detailed study which we leave to the future. 
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A Preliminary scan of the l^l~6J SM backgrounds 

In this Appendix we describe a preliminary background scan used to select the relevant back- 
grounds for the H signal discussed in Section 14.51 In this scan, we computed cross sections for 
all SM processes with QJl^l^ final state using the kinematical cuts Eq. fl4.4p as well as the final 
state invariant mass cut 

500GeV < mtot,inv < 600GeV . (A.l) 

This invariant mass range, even though more narrow than the range (14.51) used in the full analysis 
of Sections 14.51 14.61 covers the most important region near the H mass and is suitable to select 
the relevant backgrounds. The cross sections were computed either using madgraph+decay or 
ALPGEN. The processes with cross section comparable or greater than the signal are those listed in 
Table m For example, the backgrounds ZVV2j, ZtVj, hh, hV2j, hVV, ZVVV were analyzed and 
found to be negligibly small. ZVV2j, ZVAj and hAj backgrounds are also negligible, although 
this conclusion is not immediate and relies on the use of reconstruction cuts. Below we give a 
more detailed discussion of these three cases. 
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Process 


0"n=l 


0"n=2 


0"n=3 


C"n=4 


ZWnj 


2357(3) fb 


1550(3) fb 


418(1) fb 


< 200 fb 


ZZnj 


479(1) fb 


280(2) fb 


67.2(4) fb 


<30fb 



Table 5: ZVnj cross sections. The last cohimn is an extrapolation obtained taking ^ a" -I ~^ ^"ig j ^ 
a naive suppression factor for the step from n = 3 to n = 4. 

A.l Mj 

A full matrix element calculation of this process summing over gluon fusion(GF) and vector boson 
fusion(VBF) is not available. Separate computations of GF and VBF are available in madgraph, 
and GF is also available in alpgen. We used the faster alpgen for the GF case, although it has 
a caveat that all jets are assumed light and at least one has to be a gluon. Enforcing cuts Eq. 
(14.41) . setting jjj^ = + Y^jPt ' ^^"^ with no final state invariant mass restrictions we got 

a{gg ^ hg3j) = 253(1) fb. (A.2) 

The VBF was computed with madgraph by setting the ggh effective coupling to zero. For 
jj;^ = m| and with no final state invariant mass restrictions we got a cross section a factor of 5 
smaller than (lA.2p . We concluded that the GF process alone should give a reasonable estimate 



for the total cross section, and in particular that it is consistent to neglect interference effects. 

We thus proceeded to generate a sample of hgSj events using the alpgen GF process. Then, 
using DECAY, the h was made to decay into Enforcing cuts Eqs. ( 14.41) . ( 14.71) and ( lA.lll 

yields a = 0.42(2) fb, at this stage comparable with the signal (see Eq. ( 14.3p ). We then subjected 
this sample to the reconstruction cuts Rl and R2 discussed in Section The rejection efficiency 
was found to be ~ 1/30, reducing this background to a level which can be safely neglected. 

A.2 ZVV2j 

We will only discuss ZWW2j, which is the largest of the ZVV2j backgrounds. Its cross section 
after cuts Eq. (14.41) and Eq. (lA.ip was computed through ALPGEN and amounts to 4.33(2) fb. 
Decaying the Zs into leptons and Ws into quarks, we get 

a{pp 2j4:Jl+r) = 0.2 fb. 

This already looks negligible, but for reasons which will be clear in the next section we generated 
a sample and computed the rejection efficiency for the reconstruction cuts Rl and R2 of Section 
14.61 The found rejection efficiency was ~ 1/50, making this background completely negligible. 

A.3 ZVij 

Simulation of this process with a generic matrix element generator like madgraph would require 
computational resources out of our reach. At the same time we are not aware of any specialized 
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code simulating this process. For this reason we resorted to an order-of-magnitude estimate for 
this background, which seems to indicate that it is subleading. 

First of all, we simulated ZWnj and ZZnj in alpgen until n = 3 (the current limit). En- 
forcing production cuts Eq. fl4.4p and Eq. flA.ip we obtained the results given in Table Conser- 
vatively extrapolating these numbers, we got an estimate a {ZV4j) < 230 fb. Making the Z and 
W decay, we get 

a{pp ZVAj 2jUri^) < 11 fb. (A.3) 

quite a bit higher than the signal total cross section of Eq. (14.31) . However, the background is 
reducible, and we would like to estimate the rejection efficiency for the reconstruction cuts Rl 
and R2 of Section ITUl We cannot produce a sample of ZVij, but we can use the intuitively clear 
fact that ZV4j is more reducible than ZWW2j analyzed in the previous section. Thus we expect 
that the rejection efficiency should not be worse than 1/50 found for ZWW2j. From ( 1A.3I) . this 
gives a < 0.21 fb for the events which passed the reconstruction cuts. Comparing this bound 
with values reported in Table [2] we conclude that adding this channel would enhance the total 
background cross section by no more than 20%. As such, this would not change the conclusion 
about the observability of H. 
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